###########################CARREGAR PACOTES##############################
#########################################################################
require(RcmdrMisc) #pacote para carregar arquivos do Excel
require(agricolae)
require(dae)
require(ExpDes)
#Mudar diretório
setwd("D:/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
setwd("C:/Users/ander/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")

source("funcoes.R")

#########################################################################
#####################SUBSTITUIÇÃO PROTEINA###############################
#########################################################################


##Importar dados
dados=readXL("coelhos.xls",sheet = 1,stringsAsFactors=TRUE)
attach(dados)
str(dados)

trat=as.factor(sub)
####GERAR MODELO
model=aov(peso~trat)

#####PRESUPOSIÇÕES DA ANOVA
#Gerar os erros padronizados do modelo
e=rstandard(model)

#Teste de Shapiro-Wilk - NORMALIDADE
shapiro.test(e)

#Teste de Levene - Homogeneidade de Variância
leveneTest(model)

##Teste de Durbin-Watson - Independencia
set.seed(12345)
durbinWatsonTest(model)

          
####Anova com Teste de Tukey####################
crd(sub, peso, quali = FALSE,  sigT = 0.05, sigF = 0.05)

grafico.reg(sub,peso,grau=1,xlab="Substituição de Proteina",ylab="Peso")

#########################################################################
#######################################HCG###############################
#########################################################################


##Importar dados
dados=readXL("utero.xls",sheet = 1,stringsAsFactors=TRUE)
attach(dados)
str(dados)

trat=as.factor(hcg)
####GERAR MODELO
model=aov(peso~trat)

#####PRESUPOSIÇÕES DA ANOVA
#Gerar os erros padronizados do modelo
e=rstandard(model)


#Teste de Kolmogorov-Smirnov - NORMALIDADE
m=mean(e)
dp=sd(e)
ks.test(e,"pnorm",sd=dp)

#Teste de Levene - Homogeneidade de Variância
leveneTest(model)

##Teste de Durbin-Watson - Independencia
set.seed(12345)
durbinWatsonTest(model)

          
####Anova com Teste de Tukey####################
crd(hcg, peso, quali = FALSE,  sigT = 0.05, sigF = 0.05)

grafico.reg(hcg,peso,grau=2,xlab="Substituição de Proteina",ylab="Peso")